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Abstract 

We report branching fractions and charge asymmetries for exclusive decays of charged and neutral 
B mesons to two-body final states containing a charmonium meson, J/'j/j or ^(2S). This result is 
based on a 29.4 fb“^ data sample collected at the T(4S) resonance with the Belle detector at the 
KEKB asymmetric e^e“ collider. 

PACS numbers: 13.25.Hw,14.40.Gx,14.40.Nd 


*on leave from Nova Gorica Polytechnic, Nova Gorica 
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I. INTRODUCTION 


Investigation of UP violation is one of the key issues facing elementary particle physics. 
Recently, BaBar|^ and Belle[0 have observed large time-dependent UP asymmetries in the 
neutral P-meson system[^. Decay modes of neutral B mesons to hnal states containing 
charmonia were used for these measurements due to their clean experimental signatures 
and straightforward theoretical interpretation. It is expected for the same reasons that 
exclusive charmonium modes will continue to play a major role in UP studies, with rarer 
modes contributing as the body of data grows in magnitude and different aspects of the UP 
question move to the forefront. For example, the Kobayashi-Maskawa model[0 predicts small 
direct UP violation for B J/ipK^ and B Large direct UP violation would 

indicate new physics|^]. In addition, the dominant mechanism for charmonium production 
in P-meson decay is color suppressed, so precise measurements of rates to the exclusive 
modes can provide important information toward the understanding of color suppression. 

In this paper we report measurements of branching fractions and charge asymmetries for 
the exclusive decays of P mesons to the two-body hnal states xjjh, where is J/V' or '0(2S) 
and h is one of the light mesons tt^, or 7r°. We used a 29.4 fb~^ data set which 

contains 31.9 million PP events collected with the Belle detector[^ at KEKB|p. 


II. THE BELLE DETECTOR 

KEKB is an asymmetric electron-positron storage ring that collides 8.0 GeV electrons 
with 3.5 GeV positrons at the T(4S) resonance (10.58 GeV center-of-mass energy). The 
T(4S) resonance is boosted by P'y = 0.425. There is a 22 mrad crossing angle between the 
electron and positron beams at the interaction point. 

The Belle detector surrounds the beam crossing point. It is a large solid angle spec¬ 
trometer with a 1.5 T superconducting solenoid magnet. Gharged particles are detected 
by a three layer double-sided silicon vertex detector (SVD) and a 50 layer cylindrical drift 
chamber (GDG) hlled with a helium-ethane gas mixture. The tracking acceptance covers 
the laboratory polar angle between 6 = 17° and 150° {z is along the beam direction), corre¬ 
sponding approximately to 92% of the full solid angle in the center-of-mass frame (GM). The 
resolutions in impact parameter and momentum are measured to be 55 fim for a 1 GeV/c 
charged particle and CTpJpt = (0.30//3 © 0.19pt)%, where pt is the transverse momentum 
in GeV/c. A GsI(Tl) electromagnetic calorimeter (EGL) is located inside the solenoid coil 
and covers the same solid angle as the charged particle tracking system. It detects electro¬ 
magnetic showers with a resolution of cte/P = (1.3 © 0.07/P © 0.8/P^/^)%, where P is in 
GeV. 

Gharged hadron identihcation is accomplished by combining the response from an array of 
1188 silica aerogel Cerenkov counters (ACC), an array of 128 time-of-flight counters (TOE) 
and specihc ionization [dE/dx) measurement in the GDC. An iron flux-return yoke outside 
the solenoid is comprised of 14 layers of 4.7 cm-thick iron plates interleaved with a system 
of resistive plate counters (KLM) that are used for muon identihcation. The Belle detector 
is described in detail elsewhere [||. 
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III. EVENT SELECTION 


Hadronic events are selected by requiring (1) at least three reconstructed charged tracks, 
(2) a total reconstructed ECL energy in the CM in the range 0.1 to 0.8 times the total CM 
energy, (3) at least one large-angle cluster in the ECL, (4) a total visable energy (sum of 
charged tracks and neutral showers not matched to tracks) greater than 0.2 times the total 
CM energy, (5) absolute value of the 2 ; component of the CM momentum less than 50% of the 
total CM energy, and (6) a reconstructed primary vertex that is consistent with the known 
location of the interaction point. These selection criteria are determined by Monte Carlo 
simulation to be 99% efficient for signal events. To suppress two-jet non-T(4S) background 
relative to BB events we require that R 2 < 0.5, where R 2 is the ratio of the second to zeroth 
Fox-Wolfram moments [^. To remove charged particle tracks that are badly measured or 
do not come from the interaction region, we require dz < 5 cm for all tracks other than 
those identihed as decay daughters of K^, where dz is the absolute value of the coordinate 
along the beam direction at the point on the track nearest the origin. 

The decay modes considered are listed in Table |. (Hereafter the inclusion of the charge 
conjugate states is implied.) 


TABLE I: Analyzed decay chains. 


Primary mode 

Secondary mode(s) 


B~ J/ipTr~ 

J/iP l+l- 


—> J/lpTT^ 

J/lp —> /■*■/“, vrO —> yy 


B- J/ipK- 

J/lP l+l- 


50 ^ J/ipK^ 

J/lp —> Bl~,Kg 7r~^7r~ 


B- iP{2S)K- 

ip{2S) —> l~^l~, ip{2S) J/ip'K~^'K~{J/Ip - 

-4 l+l-} 

BO ^ ipi2S)K^s 

ipi/2'S,) , ip{2S) J/ipTT^Tr~{J/Ip - 

-4 l+l-}, vr+vr- 


A. J/ip{'ll; {2S)) Candidates 

In this analysis, J/'0 candidates are reconstructed from oppositely charged lepton pairs, 
or e^e“. Lepton candidates are selected with tight or loose criteria depending on 
the background level for each mode. For muon tracks, tight identihcation is based on track 
penetration depth and hit scatter in the KLM system, while loose identihcation requires 
that the track have an energy deposit in the ECL that is consistent with that of a minimum 
ionizing particle. Electron tracks are tightly identihed by a combination of dE/dx from the 
CDC, E/p {E is energy deposit in the ECL and p is momentum measured by the SVD and 
the CDC), and shower shape in the ECL. For weak electron identihcation, either dE/dx or 
E/p is required to be consistent with the electron hypothesis. 

For the identihcation oi J/'ip dilepton decays in the B —> J/ipK modes we require one 
tightly and one loosely identihed lepton. For the ip{2S)K and J/ipn modes, both lepton 
candidates are required to be tightly identihed. We correct for hnal state radiation or 
bremsstrahlung in the inner parts of the detector by including the four-momentum of every 
photon detected within 0.05 radians of the original electron or positron direction in the 
invariant mass calculation. Figure |l| shows the dilepton mass distribution near the J/ip mass. 
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The mass resolutions are 9.3 MeV/c^ and 10.6 MeV/c^ in the peak region for and 

respectively. Since there are still small radiative tails, as can be seen in Figure 0, we 
use asymmetric invariant mass requirements, (—60 < m^+^- — mj/^(^^( 2 S)) < 36) MeV/c^ 
and (—150 < rue+e- “ '^j/y(y( 2 S)) < 36) MeV/c^, for the and e’''e“ pairs respectively. 



FIG. 1; The invariant mass distributions for (a) and (b)e'''e . In these figures, both leptons 

are tightly identified. 

To identify '0(2S) —> J/'i/ivr+7r“ candidates, we combine J/il) candidates with pairs of 
oppositely charged tracks that have a tt+tt” invariant mass greater than 400 MeV/c^. The 
'0(28) and J/'ip candidates’ mass difference is required to be consistent with the known 
difference, (0.58 < - 07 .;+;-^+,^- — rrii+i- < 0.60) GeV/c^. This range corresponds to ±3(T in 
detector resolution. Figure ^ shows (a) the invariant mass distribution of '0(28) —>■ ZZ, and 
(b) the mass difference of — mu. 


B. Light meson candidates 

In the analysis for B~ —»• J/'0('0(28))iF“, all charged tracks (other than those used for 0 
reconstruction) are used as kaon candidates in order to eliminate the systematic error from 
particle identihcation. This does not introduce any serious additional background, because 
the principal background is expected to be from B~ —>• J/0(0(28))7r“ decays, which occur at 
a much lower rate than J/'0('0(28))iF“. The prompt charged pion candidates are conversely 
required to be strongly identihed as pions {P{tt/K) > 0.9), where the likelihood ratio for a 
particle to be a charged pion, P{7i/K) = Proh^Ti)/{Proh^n) + Prob{K)), is calculated using 
dE/dx measured in the GDC and the response of the AGG. 

For the analysis of neutral B meson decays, the reconstruction of —> tt+tt" is made 

by selecting pairs of oppositely charged tracks with tt+tt" invariant mass between 482 and 
514 MeV/c^. This criterion retains 99.7% of —>■ vr+vr” decays with detected tracks, 

based on a double Gaussian £t to the mass peak of the data (the average mass resolution is 
4.4 MeV/c^.). In order to reduce combinatorial background further, we require that: 
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FIG. 2: (a) The invariant mass distribution of 'ip{2S) —> candidates, (b) the mass difference 

of mi+i-^+^- — rrii+i-. 


• if both pions have associated S VD hits, the points of nearest approach of the two tracks 
in the projection onto the plane perpendicular to the beam line (r- 0 ) are separated in 
the beam direction (z) by less than 1 cm, 

• if only one of the two pions has associated SVD hits, the distance of nearest approach 
to the interaction point in the r-0 projection be greater than 250 pm for both tracks, 

• if neither of the two pions have associated SVD hits, the 0 coordinate of the 7 r’'' 7 r“ 
vertex point and the 0 direction of the tt+tt" candidate’s three-momentum agree within 
0.1 radians. 

The Kg hnding efficiency after track selection is 95%. 

In the selection of > J/ 07 r°, the high momentum vr^’s are reconstructed from pairs of 
detected photons. The invariant mass is required to be 118 MeV/c^ < < 150 MeV/c^ 

(mass resolution is 5.3 MeV/c^). The 7 r° candidate is also to have a good mass constrained 
£t. 


C. B Meson Reconstruction 

B mesons are reconstructed by combining a charmonium meson candidate with a kaon 
or pion candidate, as described above. The energy difference, AE = i^cand — -^beam, and 
the beam-energy constrained mass, Mbc = \J— P^and) are used to separate signal from 
background (i^beam is the beam energy, T^cand and Pcand are the B candidate energy and 
momentum, all calculated in the T(4S) center of mass frame). 

In this calculation, kinematic fits are performed with (1) mass and vertex constraints 
for the J/0 or '0(28) di-lepton decays and decays, and (2) a mass constraint for the 
0(28) —> J/07r+7r“ and 7r° —>• 77 decays. Figure ^ shows the distribution for B 
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5.200 5.225 5.250 5.275 5.300 -0.20 -0.10 0.00 0.10 0.20 

Beam Constrained Mass (Gev/c^) AE (GeV) 


FIG. 3: The distribution of (a) AE versus M^c and (b) AE for B —> The background 

from B —> is seen at lower AE, while that from B —> J/'i/ivr^ is at higher AE. 


candidates in the M\yc-AE plane as well as in AE after projecting ont candidates with Mbc 
between 5.27 and 5.29 GeV/c^. In the hrst plot an excess of candidates is clearly apparent 
in the signal region, indicated by the rectangle. 

In order to determine yields, we £t the Mbc distribntions after applying the following 
reqnirements on AE\ for all modes except J/V’vr, (—40 < AE < 40) MeV; for the J/'ip'K^ 
decay mode, (—100 < AE < 50) MeV, as the AE distribntion has a long tail at negative 
valnes dne to material in the detector and energy leakage; for the J/'07r“ mode, (—10 < 
AE < 40) MeV, to snppress a backgronnd from B —>• J/ipK dne to misidentihcation of 
K~ as 7r“. 

The £t of the Mbc distribntion is performed with the snm of a Ganssian for signal and 
the ARGUS fnnction|[Tl[] for backgronnd (Fignre The resolntion in Mbc is dominated by 
the energy spread of KEKB. We test the resolntion agreement between MG and data nsing 
the mode B~ —»• J/'ipK~ (which has the highest statistics). The agreement is very good 
and we nse the MG predicted widths for each mode to £t the Mbc histograms. The signal 
yield and the normalization of the backgronnd are allowed to vary in the £t. The results are 
shown in Table |^. 

For the J/V’tt” mode, as shown in Figure the background from B~ J/ipK~ peaks in 
the signal region of Mbc but accumulates near AE ~ —70 MeV due to kinematic differences 
from the signal mode. To insure that it does not intrude into the signal in the Mbc ht, a £t is 
performed on the AE distribution with two separated Gaussians and a hrst-order Ghebyshev 
polynomial function (Figure ^). The signal yield obtained from the AE distribution fit is 
consistent with the Mbc fit- The background yield from J/ipK~ is also consistent with the 
expectation from the mis-identification probability. 
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Beam Constrained Mass (Gev/c^) 


FIG. 4: The distribution of Mbc for (a) B —> , (b) J/tl^Kg, (c) B 

'ip{2S)K-{il:{2S) l+l-}, (d) B^ V’(2S)i^g{V’(2S) ^ /+/“}, (e) B- ^jJ{2S)K-{^P{2S) 

J/'t/tyr+TT”}, (f) B^ '^jJ{2S)Kg{^p{2S) — > Jf'ip'jr~^TT~}, (g) B~ — t- J/'ijjTT~ and (h) B^ — > 

IV. RESULTS 

A. Branching Fractions 


The reconstruction efficiencies are determined by Monte Carlo simulations (MC) based 
on GEANT0 and are listed in Table The number of BB events is measured to be 
(31.9 ± 0.3) X10®. In the calculation of the branching fraction, the production rates of 
B^B~ and B'^B^ pairs are assumed to be equal. We use the secondary branching fractions 
listed in Table |T1| [O. The resulting branching fractions for each reconstructed decay 


chain are summarized in Table |1|, where the first errors are statistical and the second are 
systematic. The measurement values for the two '0(2S) modes of i? —^ 'ip{2S)K are consistent 
within their errors and the combined results are also listed in the table (taking into account 
correlated and uncorrelated errors). 

The sources of systematic error are shown in Table Eg The dominant uncertainty arise 
from the uncertainty in the tracking efficiency. 

The tracking efficiency uncertainty is determined to be 2% per track from a comparison 
of the yields for rj —>■ Tr+vr^Tr® and r; —>■ 77 with MC expectations. The pion tracks from 
decay are from a displaced vertex and thus may have a larger systematic error. We include 
a 3.5% per track uncertainty for these tracks (see below). 

The uncertainty in the selection efficiency is checked by comparing yields for a sample 
of high momentum —>■ tt+tt” decays before and after applying the selection criteria. 

The efficiency difference between data and the Monte Carlo simulation is less than 1.0%. 
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AE (GeV) 


FIG. 5: The /S.E distribution for —>■ The signal peak is seen around zero. The peak at 

—0.07 GeV/c^ is from —> J/'ipK^. In this figure, we require 5.27 < Mbc < 5.29 GeV/c^ 


As one cross-check of reconstruction, we also estimate the difference between data 
and non-T(4S) MC directly, using the yield ratio between and K~n^ 

with from D* D^tt decay. In this case, with momentum higher than 3.0 GeV/c 
are selected. The difference of the ratio between data and MC is also smaller than 1%, where 
a large systematic error arises from the uncertainties of the world averages for the branching 
fractions. 

The high momentum efficiency is checked by taking the ratio between K^7r~7i^ 

and K^n~ with high momentum D^’s. D^’s generated from D* decay with a slow 

pion are selected. We assign a 7% uncertainty to the 7r° efficiency. 

The efficiency of lepton identihcation is checked by comparing the J/ijj yield with one 
lepton tightly identihed against the yield where both leptons are tightly identihed. We hnd 
that the efficiencies for tightly identihed electrons and muons are 96% and 94%, respec¬ 
tively. The systematic errors from lepton identihcation are determined to be 2% per tightly 
identihed lepton. The error for loosely identihed leptons in negligable. 

For the B~ mode, the identihcation of high momentum charged pions is studied 

by comparing D*^ , where K~7r~^, between data and MC. In this decay mode, 

the mass peak is reconstructed with small background without any particle identihcation 
requirements. The systematic uncertainty is determined by examining the diherence in yield 
before and after applying particle identihcation. We assign a systematic uncertainty of 2% 
to the pion identihcation efficiency. 

We also study the systematic error arising from the background in the ht of the Mbc 
distribution. The ARGUS function represents the Mbc distribution for the AE sidebands 
well. However there may be background decay modes that peak in the signal region. We 
checked this with inclusive J/'ip MC and hnd no evidence for peaking background. 
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TABLE II: Signal yields and branching fractions for each mode. Signal yields are determined by fit¬ 
ting. The errors are statistical (first error) and systematic (second), except for the combined 'ijj{2S) 
modes where the total error is listed. Efficiencies for modes with mesons are for reconstructing 
B 

Decay Mode Yield Efficiency(%) B. E. (xlO“^) 

43.9 ± 6.8 33.3 0.38 ± 0.06 ± 0.03 

24.0 ± 5.0 27.2 0.23 ± 0.05 ± 0.02 


B-^J/ipK- 2102 ±46 55.3 10.1 ±0.2 ±0.7 

5° ^ J/ipK^ 453 ± 21 30.5 7.9 ± 0.4 ± 0.9 


B- 

^p{2S)K- 




6.9 ±0.6 



V’(2S) ^ l+l- 

173 ± 13 

51.6 

7.3 ±0.6 ±0.7 



-0(28) —> J/f/iTr+TT” 

170 ± 13 

23.2 

6.4 ±0.5 ±0.8 

BO ^ 

V^(2S)iL0 




6.7 ± 1.1 



V^(2S) ^ l+l- 

38.5 ± 6.2 

27.5 

6.1 ±1.0 ±0.8 



'i/i(2S) ^ J/iIjtt+'k- 

51.2 ±7.2 

12.0 

7.4 ± 1.0 ± 1.3 

TABLE III: Branching fractions used for secondary charmonium decay s|[l^]. 

Decay mode 




Branching fraction 

J/V-- 

-> e^e 




0.0593 ± 0.0010 

J/V-- 

fi+iJ,- 




0.0588 ± 0.0010 

V^(2S) 





0.0073 ± 0.0004 

^(2S) 





0.0070 ± 0.0009 

V^(2S) 

i 

+ 

1 




0.305 ± 0.016 


B —> J/'lfjTT 
B^ —> J/'07r^ 


B. Charge Asymmetries 

The yields for positive and negative B mesons decays are measnred separately nsing the 
method described above. The charge asymmetries, dehned by 

^ Br{B~ —> charmonium ± K~{'k~)) — Br{B^ —> charmonium ± 

Br{B~ —> charmonium ± ± Br{B^ charmonium ± iL+(7r+)) ’ 

(see [^) are calculated assuming the same efficiencies for both charged decays. The results 
are shown in Table 0. The efficiency difference between positive and negative particles 
is determined by using 3.96 x 10^ and 3.33 x 10^ events for and 

K~7i^/D^ —> K~^7i~ decays, respectively. We calculate the efficiency ratios /€.„-+ = 
1.011 ± 0.015 and ^k- = 1-004 ± 0.017 using the following formulas: 

e^- _ N{D- K+'k-t:-)N{D^ K-'k+) 

7^ ~ N{D+ K-'k+t:+)N(W AT+vr-) ’ ^ ^ 
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TABLE IV: The dominant sources of systematic errors (in %). 


Decay mode 

Tracking 

Lepton Hadron 

Kg{7r^) Charmonium 

Monte 

Total 


Efficiency 

ID ID 

Efficiency Branching 

Carlo 




Efficiency Efficiency 

Fractions 

Stats. 


B~ —> J/'07r“ 

6.0 

4.0 2.0 


1.2 

1.7 

7.8 


4.0 

4.0 

7.0 

1.2 

1.9 

9.3 

B- J/ipK- 

6.0 

2.0 


1.2 

1.4 

6.6 

B^ 

11.0 

2.0 

1.0 

1.2 

1.8 

11.4 

B- ^^{2S)K- 







V’(2S) ^ l+l- 

6.0 

4.0 


5.0 

1.4 

8.9 

V’(2S) —>■ 

10.0 

4.0 4.0 


5.3 

2.1 

12.8 








V’(2S) ^ l+l- 

11.0 

4.0 

1.0 

5.0 

1.9 

12.9 

'0(28) —> J/lpTT+TT- 

15.0 

4.0 4.0 

1.0 

5.3 

2.9 

17.1 

TABLE V: Charge asymmetry for each mode. Errors are statistical only. 

Decay mode 


Yield(-) 

Yield(-I-) 




B+ —> J/07r^ 


21 ±5 

22 ±5 


-0.023 ±0.164 

B+ J/'iIjK+ 


1024 ± 32 

1078 ± 33 


-0.026 ± 0.022 

B+ 0(2S)(/+/-)A± 


79 ±9 

93± 10 


-0.081 ± 0.078 

B+ 'il>{2S){J/'tp7T+7r- 

-)K+ 

68 ±8 

102 ± 10 


-0.200 ± 0.075 

Total iB+ J/0(0(2S))A:±) 

1171 ± 34 

1273 ± 36 


-0.042 ± 0.020 


eK- N{D- K+7r-7r-)N{D^ K-ir+f 

7 ^ ~ N{D+ K-'k+t:+)N(W K+Ti-y ^ ’ 

No significant efficiency differences are observed for either pion or kaon tracks. Thus, we 
do not correct the central values but we do include the error of the efficiency differences in 
the systematic errors. 

Finally, we hnd the charge asymmetries —0.023±0.164±0.015 and —0.042±0.020±0.017 
for the charmonium+TT mode and the charmonium+iT mode, respectively. Our results are 
consistent with zero asymmetry and previous measurements p!4| . 


V. CONCLUSION 

We have reported measurement of B meson branching fractions to two-body hnal states 
that include a J/ijj or "0(28) meson and a 7r° or vr^. A total of 31.9 million BB 

events accumulated at the T(4S) resonance are used for this analysis. Our results are in 
good agreement with previous measurements |T5|, pffif . Charge asymmetries are also measured 
and found to be consistent with zero. 
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